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1. Introduction 

Overlap fermions provide the opportunity to investigate QCD at small quark masses. Be¬ 
ing computationally extremely expensive, no dynamical results on physically relevant lattices are 
available yet. In this work, we present spectrum and quark mass results of a quenched study. 
We calculate light and strange quark masses, and we check several predictions of quenched chiral 
perturbation theory. Results for nucleon matrix elements are presented in [[I]]. 

2. Action, lattices and the scale 

Our overlap operator is 

D = ( l -^f) DN + mq ’ ^ = ^(l + X(Xtx)- 1/2 ), (2.1) 

with the Wilson kernel operator X = D\y — £ We use a polynomial approximation (see [2j]) for 
X (X^X) , the degree of which is adjusted properly to get an overall residual in the inversion 

of the overlap better than 10 7 . The value of p is chosen to satisfy the requirements of low 
condition number and good localisation properties. A way to improve the condition number is the 
exact treatment of the lowest eigenvalues. Our choice, p = 1.4, is a compromise between small 
condition number and large localisation mass (i.e. we aim at the upper left hand corner in Fig. [I]). 
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Figure 1: Condition number and localisation 
mass (1-norm) for /3 = 8.45, 10 eigenvalues pro¬ 
jected out. 



Figure 2: Quenched chiral fits to f n . 


We use the 1-loop tadpole improved Luscher-Weisz gauge action [|J]. 


S[U) = 


co £ ^ReTr (1 - U v ,i aq ) +d £ ^ReTr (1 - U Kct ) + c 2 £ ^ReTr(l - t/ par ) 

plaq ^ rect ^ par ^ 


( 2 . 2 ) 


The lattice gauge coupling f5 is defined as /3 = 6co/g 2 , and the values of c\ and C 2 are fixed by the 
value of co and the plaquette expectation value uq. The benefits of the Luscher-Weisz compared 
to the Wilson gauge action are a better condition number of (as long as no gauge smearing is 
involved) and absence of dislocations. 
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We use Jacobi smeared sources (jc = 0.21, n = 50). To set the scale we compute 


fn{m q ) 


m q A s i 



(2.3) 


for each of our quark masses, where m n and A s i jSS are the mass and amplitude from cosh fits to the 
smeared-local and smeared-smeared pseudoscalar correlators. We extrapolate the lattice numbers 
to the chiral limit by the leading order formula 


fnim q ) =f K 


1+4 L 


q m l 

'5 f 2 
J 71 


(2.4) 


The resulting values of a are given in Table [I] (for the values of L| see the legend of Fig. |2|). Note, 
that the lattice constant is about 10% larger than in 0 where the Sommer scale was used. Such a 
scale ambiguity is a familiar phenomenon of quenched simulations. 


j8 

lattice 

cfg. 

a/ fm 

a -1 / MeV 

8.0 

16 3 x 32 

230 

0.153(3) 

1290(30) 

8.45 

16 3 x 32 

250 

0.105(2) 

1870(40) 

8.45 

24 3 x 48 

200 

0.105(1) 

1870(20) 


Table 1: Lattices used in the simulation. The scale was determined via f K (see Fig. |j). 


3. Quark masses 

The best signal for the pseudoscalar mass is obtained from the correlator of the time compo¬ 
nent of the axial current. An example of the results is displayed in Fig. [3]. The quenched chiral 




Figure 3: The squared pseudoscalar mass as Figure 4: A plot of m\/m q reveals the existence 

of quenched chiral logs. 


function of m q together wit the fit Eq. (3.2). 
perturbation theory prediction for m K reads 

i 2 n =Am q ^1 — 5 ^ln ^Am q /A^j + lj^ + G (m q ) 


m. v 


(3-1) 
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lattice 

8 

8.00, 16 3 x 32 

0.16(4) 

8.45, 16 3 x 32 

0.3(1) 

8.45, 24 3 x 48 

0.15(3) 


Table 2: The chiral log parameter 8. 


where the quenched chiral log appears with a prefactor 8 giving rise to a singularity in m 2 n /m q . The 
fit is shown in Fig. [|. Fixing A x to 1 GeV, we obtain the values of 8 in Table |2| Although the chiral 
logs will strongly affect the determination of ni£ = 1 / 2 (m u + mf), we will ignore this problem in 
the following. 

The fit functions we use to determine the bare light and strange quark masses are 


9 9 

m n = Ante + Bni£ In me +Cni(, 


m~ K = A 


ni( + m s n m £ + m s / m s In m s — me In me 


+ B 


-1 


m s — me 


+c ( 


me + m s 


(3.2) 


The renormalisation is done non-pertubatively 1 in the Rl' —MOM scheme with a variant of the 
method introduced in 0, using momentum sources |7|]. This variant requires only small statistics, 
and one does not need to perform inversions for each operator under investigation (although one 
needs to invert the Dirac operator for each momentum separately). We implement the renormalisa¬ 
tion condition for the vertex 

Tr(rg n r o 1 Born )=12, (3.3) 

on the lattice. Putting in the definition of the renormalised vertex we obtain 


= Zs = 


ikTr(r s ry^ ra ) 


(3.4) 


In the following, we use Aq for the denominator of the r.h.s of Eq. (3.4) for O = {A. V. P..S}. The 
wave function renormalisation Z q is computed from Z q = Z 4 A 4 , where we use the axial Ward 
identity to compute Z ,\, so that it can be extrapolated 


Za 


2m AP (t)P( 0)) 

lim lim 

(5 4 A 4 (t)P (0)) 


(3.5) 


The chiral extrapolation is illustrated in Fig. |5| 

We now proceed with the computation of A^and As. We exploit the identity is A 5 = Ap, a 
consequence of chiral symmetry. On the lattice both are modified by zero mode effects oe 1 /( am q ) 2 , 
which are finite volume artefacts. Additionally, spontaneous chiral symmetry breaking produces a 
1 /( 1 am q ) contribution to Ap tt . Thus, we extract As as the common constant from fits 


Ap" = -^ + -^-+A s + PAcrm 2 
a L niq am q 1 

A^" = — 2—2 sja 2 m q . (3.6) 

Am- H 

1 A perturbative calculation of the renormalisation factors can be found in |^J. 
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Figure 5: Z A from the axial Ward identity. Figure 6: Z RGI for both values of jS 


for each value of jl 2 . A a is computed in a similar fashion from linear fits to A |[ 11 and Ay tt with a 
common constant. In the next step we need to identify the region where non-perturbative effects 
are small. To this end we remove the renormalisation group running from the RI' — MOM renor¬ 
malisation constant Z$ = Z a Aa/A$ using the 4-loop results from |[s]|. The result is shown in Fig. g. 
where the linear regions are clearly identified. We extract the intercept and quote the results in the 
MS scheme at f = 2 GeV in Table g. Note, that the lattice artefacts (represented by the slope) 


p 

zf G/ 

slope 

Z(p( 2 GeV) 

8.00 

1.18(2) 

-0.007(2) 

0.85(2) 

8.45 

1 . 02 ( 1 ) 

-0.004(1) 

0.73(1) 


Table 3: Z s 


P 

mf 

mp 

8.00 

4.2(1) 

127(1) 

8.45,16 3 32 

3.5(2) 

118(3) 

8.45,24 3 48 

4.1(1) 

119(1) 


Table 4: Light and strange quark masses 


scale with a 2 , as expected for (the automatically A (a) improved) overlap fermions. The resulting 
quark masses are given in Table Q The light quark mass scales very well, but we see a volume 
dependence. For the strange quark mass we see some dependence on a, but not on the volume. 


4. Vector meson and nucleon masses 

In this section we extract vector meson and nucleon masses from our data and compare their 
chiral behaviour to predictions of quenched chiral perturbation theory. The vector meson mass 
(Fig. 0) confirms the prediction of a negative value of Cy/2 in a fit according to m p {m n ) = m p + 
Ci/2 m jr + C\m 2 K + The prediction for the nucleon has the same structure, but it is harder 

to confirm the negative sign of Cy/ 2 for my, since the errors make the results compatible with zero 
unless one uses the low mode averaging technique. In this case we find a value of Cy/ 2 = —0.32(18) 
at /J = 8.45, 24 3 48, which is in good agreement with the chiral perturbation theory prediction 
Cy/2 = 3/2 ti(D - 3F) 2 8 = — 0.25(5) 2 (Fig. g). 

5. Summary 

In a quenched overlap simulation we have computed the light and strange quark masses 

2 D and F are axial current matrix elements, D=0.81(3), F=0.47(4) [^J 
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m^/MeV 2 


Figure 7: Quenched chiral fits to the vector me- Figure 8: The same for the nucleon. The 
son. dashed line is a fit to low mode averaged data. 


mf s {2 GeV) = 4.1(1) MeV and (2 GeV) = 119(1) MeV. We demonstrated the existence 
of quenched chiral logs and the expected chiral behaviour of f n , m p , and my. 
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